Simulation of non-Markovian Processes in BlenX

BlenX is a programming language explicitly designed for modeling biological processes inspired by Beta-binders. The actual framework assumes biochemical interactions being exponentially distributed, i.e., an underlying Markov process is associated with BlenX programs. In this paper we relax this condition by providing formal tools for managing non-Markovian processes within BlenX

Erratum to 'Core Biopsies from Prostate Cancer Patients in Active Surveillance Protocols Harbor PTEN and MYC Alterations'[European Urology Oncology 2 (2019) 277-285].

n/

Analyzing various models of Circadian Clock and Cell Cycle coupling

The daily rhythm can influence the proliferation rate of many cell types. In the mammalian system the transcription of the cell cycle regulatory protein Wee1 is controlled by the circadian clock. Zamborszky et al. (2007) present a computational model coupling the cell cycle and circadian rhythm, showing that this coupling can lead to multimodal cell cycle time distributions. Biological data points to additional couplings, including a link back from the cell cycle to the circadian clock. Proper modelling of this coupling requires a more detailed description of both parts of the model. Hence, we aim at further extending and analysing earlier models using a combination of modelling techniques and computer software, including CoSBI lab, BIOCHAM, and GINsim

Unraveling the clonal hierarchy of somatic genomic aberrations

Defining the chronology of molecular alterations may identify milestones in carcinogenesis. To unravel the temporal evolution of aberrations from clinical tumors, we developed CLONET, which upon estimation of tumor admixture and ploidy infers the clonal hierarchy of genomic aberrations. Comparative analysis across 100 sequenced genomes from prostate, melanoma, and lung cancers established diverse evolutionary hierarchies, demonstrating the early disruption of tumor-specific pathways. The analyses highlight the diversity of clonal evolution within and across tumor types that might be informative for risk stratification and patient selection for targeted therapies. CLONET addresses heterogeneous clinical samples seen in the setting of precision medicine. Electronic supplementary material The online version of this article (doi:10.1186/s13059-014-0439-6) contains supplementary material, which is available to authorized users

ETS-related gene (ERG) undermines genome stability in mouse prostate progenitors via Gsk3β dependent Nkx3.1 degradation.

21q22.2-3 deletion is the most common copy number alteration in prostate cancer (PCa). The genomic rearrangement results in the androgen-dependent de novo expression of ETS-related gene (ERG) in prostate cancer cells, a condition promoting tumor progression to advanced stages of the disease. Interestingly, ERG expression characterizes 5-30% of tumor precursor lesions - High Grade Prostatic Intraepithelial Neoplasia (HGPIN) - where its role remains unclear. Here, by combining organoids technology with Click-chemistry coupled Mass Spectrometry, we demonstrate a prominent role of ERG in remodeling the protein secretome of prostate progenitors. Functionally, by lowering autocrine Wnt-4 signaling, ERG represses canonical Wnt pathway in prostate progenitors, and, in turn, promotes the accumulation of DNA double strand breaks via Gsk3β-dependent degradation of the tumor suppressor Nkx3.1. On the other hand, by shaping extracellular paracrine signals, ERG strengthens the pro-oxidative transcriptional signature of inflammatory macrophages, which we demonstrate to infiltrate pre-malignant ERG positive prostate lesions. These findings highlight previously unrecognized functions of ERG in undermining adult prostate progenitor niche through cell autonomous and non-autonomous mechanisms. Overall, by supporting the survival and proliferation of prostate progenitors in the absence of growth stimuli and promoting the accumulation of DNA damage through destabilization of Nkx3.1, ERG could orchestrate the prelude to neoplastic transformation

Allele-informed copy number evaluation of plasma DNA samples from metastatic prostate cancer patients: the PCF_SELECT consortium assay.

Sequencing of cell-free DNA (cfDNA) in cancer patients' plasma offers a minimally-invasive solution to detect tumor cell genomic alterations to aid real-time clinical decision-making. The reliability of copy number detection decreases at lower cfDNA tumor fractions, limiting utility at earlier stages of the disease. To test a novel strategy for detection of allelic imbalance, we developed a prostate cancer bespoke assay, PCF_SELECT, that includes an innovative sequencing panel covering ∼25 000 high minor allele frequency SNPs and tailored analytical solutions to enable allele-informed evaluation. First, we assessed it on plasma samples from 50 advanced prostate cancer patients. We then confirmed improved detection of genomic alterations in samples with <10% tumor fractions when compared against an independent assay. Finally, we applied PCF_SELECT to serial plasma samples intensively collected from three patients previously characterized as harboring alterations involving DNA repair genes and consequently offered PARP inhibition. We identified more extensive pan-genome allelic imbalance than previously recognized in prostate cancer. We confirmed high sensitivity detection of BRCA2 allelic imbalance with decreasing tumor fractions resultant from treatment and identified complex ATM genomic states that may be incongruent with protein losses. Overall, we present a framework for sensitive detection of allele-specific copy number changes in cfDNA

Allele-informed copy number evaluation of plasma DNA samples from metastatic prostate cancer patients: the PCF_SELECT consortium assay

Sequencing of cell-free DNA (cfDNA) in cancer patients' plasma offers a minimally-invasive solution to detect tumor cell genomic alterations to aid real-time clinical decision-making. The reliability of copy number detection decreases at lower cfDNA tumor fractions, limiting utility at earlier stages of the disease. To test a novel strategy for detection of allelic imbalance, we developed a prostate cancer bespoke assay, PCF_SELECT, that includes an innovative sequencing panel covering ∼25 000 high minor allele frequency SNPs and tailored analytical solutions to enable allele-informed evaluation. First, we assessed it on plasma samples from 50 advanced prostate cancer patients. We then confirmed improved detection of genomic alterations in samples with <10% tumor fractions when compared against an independent assay. Finally, we applied PCF_SELECT to serial plasma samples intensively collected from three patients previously characterized as harboring alterations involving DNA repair genes and consequently offered PARP inhibition. We identified more extensive pan-genome allelic imbalance than previously recognized in prostate cancer. We confirmed high sensitivity detection of BRCA2 allelic imbalance with decreasing tumor fractions resultant from treatment and identified complex ATM genomic states that may be incongruent with protein losses. Overall, we present a framework for sensitive detection of allele-specific copy number changes in cfDNA

Transcriptional Regulation Is a Major Controller of Cell Cycle Transition Dynamics

DNA replication, mitosis and mitotic exit are critical transitions of the cell cycle which normally occur only once per cycle. A universal control mechanism was proposed for the regulation of mitotic entry in which Cdk helps its own activation through two positive feedback loops. Recent discoveries in various organisms showed the importance of positive feedbacks in other transitions as well. Here we investigate if a universal control system with transcriptional regulation(s) and post-translational positive feedback(s) can be proposed for the regulation of all cell cycle transitions. Through computational modeling, we analyze the transition dynamics in all possible combinations of transcriptional and post-translational regulations. We find that some combinations lead to ‘sloppy’ transitions, while others give very precise control. The periodic transcriptional regulation through the activator or the inhibitor leads to radically different dynamics. Experimental evidence shows that in cell cycle transitions of organisms investigated for cell cycle dependent periodic transcription, only the inhibitor OR the activator is under cyclic control and never both of them. Based on these observations, we propose two transcriptional control modes of cell cycle regulation that either STOP or let the cycle GO in case of a transcriptional failure. We discuss the biological relevance of such differences

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers

On the Computational Power of BlenX

We present some decidability and undecidability results for subsets of the BlenX Language, a process-calculi-based programming language developed for modelling biological processes. We show that for a core subset of the language (which considers only communications primitives) termination is decidable. Moreover, we prove that by adding either global priorities or events to this core language, we obtain Turing equivalent languages. The proof is through encodings of Random Access Machines (RAMs), a well known Turing equivalent formalism, into our subsets of BlenX. All the encodings are shown to be correct